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Lattice dynamics and correlated atomic motion from the atomic pair distribution 

function 
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The mean-square relative displacements (MSRD) of atomic pair motions in crystals are studied as 
a function of pair distance and temperature using the atomic pair distribution function (PDF) . The 
effects of the lattice vibrations on the PDF peak widths are modelled using both a multi-parameter 
Born von-Karman (BvK) force model and a single-parameter Debye model. These results are com- 
pared to experimentally determined PDFs. We find that the near-neighbor atomic motions are 
strongly correlated, and that the extent of this correlation depends both on the interatomic inter- 
actions and crystal structure. These results suggest that proper account of the lattice vibrational 
effects on the PDF peak width is important in extracting information on static disorder in a disor- 
dered system such as an alloy. Good agreement is obtained between the BvK model calculations of 
PDF peak widths and the experimentally determined peak widths. The Debye model successfully 
explains the average, though not detailed, natures of the MSRD of atomic pair motion with just 
one parameter. Also the temperature dependence of the Debye model largely agrees with the BvK 
model predictions. Therefore, the Debye model provides a simple description of the effects of lattice 
vibrations on the PDF peak widths. 

PACS numbers: 63.20.-e, 61.12.-q, 61.10.-i 



I. INTRODUCTION 

The pair distribution function (PDF) obtained from 
the powder x-ray and neutron diffraction experiments has 
been shown to be of great valuej-.in determining the lo- 
cal atomic structure of materials.EJ The PDF results from 
a Fourier transform of the powder diffraction spcd-rum 
(Bragg peaks -I- diffuse scattering) into real-spaceJa For 
well ordered crystals, apart from technical details, this 
is similar to fitting the Bragg peaks -I- thermal diffuse 
scattering (TDS) in the powder pattern in a manner first 
discussed by Warren.H A PDF spectrum consists of a se- 
ries of peaks, the positions of which give the distances of 
atom pairs in real space. The ideal width of these peaks 
(aside from problems of experimental resolution) is due 
both to relative thermal atomic motion and to static dis- 
order. Thus, an investigation of the effects of lattice vi- 
brations on PDF peak widths is important for at least 
two reasons: first, to establish the degree to which infor- 
mation on phonons (and the interatomic potential) can 
be obtained from powder diffraction data, and, second, 
to account for correlation effects in order to properly ex- 
tract information on static disorder in a disordered sys- 
tem such as an alloy. 

In general, powder diffraction is not considered a favor- 
able approach for extracting information about phonons 
since, not only is energy information lost in the mea- 
surement, but also the diffuse scattering is isotropically 
averaged. The lattice vibrations are best described from 
the phonon dispersion curves determined using inelas- 
tic neutron scattering and high-energy-resolution inelas- 



tic x-ray scattering on single crystals. qEI Nevertheless, 
with the advent of high-energy synchrotron x-ray and 
pulsed-neutron sources and fast computers, it is possible 
to measure data with unprecedented statistics and accu- 
racy. The PDF approach has been shown to yield limited 
information about lattice vibrations in powders,B though 
the extent of which thia. information can be extracted 
remains controversial .ffltl'EJ 



Measuring powders has the benefit that the experi- 
ments are straightforward and do not require single crys- 
tals. It is thus of great interest to characterize the degree 
to which lattice vibrations are reflected in the PDF using 
simple models, such as the Debye model, in situations 
where detailed interatomic potential information is not 
available. In this paper we explore these issues by com- 
paring both measured PDFs and those calculated from 
realistic potential models with PDFs obtained through a 
single-parameter Debye model. This comparison is car- 
ried out as a function of atomic pair separation, temper- 
ature and direction in the lattice. We find that a single 
parameter Debye model explains much of the observed 
lattice vibrational effects on PDF peak widths, including 
the temperature dependence, in crystals like Ni, Ce, and 
GaAs. However, small but non-negligible deviations from 
the Debye model calculation are evident in crystal which 
needs a long-range interaction to explain anomalies in 
the dispersion curves. 
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II. CORRELATED ATOMIC MOTION IN 
REAL-SPACE 

The existence of interatomic forces in crystals results 
in the motion of atoms being correlated. This is usu- 
ally treated theoretically by transforming the problem to 
normal coordinates, resulting in normal modes (phonons) 
that arc non-interacting, thus making the problem math- 
ematically tractible. Projecting the phonons back into 
real-space coordinates yields a picture of the dynamic 
correlations. This situation can be understood intuitively 
in the following way. Figure |^ shows a schematic di- 
agram of atomic motion in three different interatomic 
force systems, each with its corresponding ideal PDF 
spectrum. In a rigid-body system. Fig. |l|(a), the inter- 
atomic force is extremely strong and all atoms move in 
phase. In this case, the peaks in the PDF are delta- 
functions. At the opposite extreme the atoms are non- 
interacting (the Einstein model) and move independently 
as shown in Fig. 0(b) . This type of atomic motion results 
in broad PDF peaks whose widths are given by the root 
mean-square displacement amplitude (v\t?}). In real 
materials, the interatomic forces depend on atomic pair 
distances, i. e., they are strong for nearest-neighbor in- 
teractions and get weaker as the atomic pair distances 
increase. In fact, these interactions are often quite well 
described with just nearest-neighbor or first- and second- 
nearest-neighbor coupling. The case of nearest-neighbor 
interactions is shown in Fig. |^(c). In this (Debye) model 
a single parameter corresponding to the spring constant 
of the nearest-neighbor interaction is used. Here, near- 
neighbor atoms tend to move in phase with each other, 
while far-neighbors move more independently. As a re- 
sult, the near-neighbor PDF peaks are sharper than those 
of far-neighbor pairs. This behavior jasias-first analyzed by 
Kaplow et al. in a series of papersElloEj for a number 
of elemental metals. 
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FIG. 1: Schematic diagram showing an instantaneous snap- 
shot of atomic positions in (a) rigid-body model (b) Einstein 
model (c) Debye model. In (a) and (b) all PDF peaks have 
the same width independent of atom separation. In (c) the 
PDF peak width increases up to the root mean-square dis- 
placement as the atomic separation increases, a is a spring 
constant. 

PDFgetxjll respectively. The experimental PDF, G(r), 
was obtained by taking the Fourier transform of S{Q): 

G{r) = f /jf ""^ Q[S(Q) - l]sinQr dQ, where Q^ax is the 
maximum momentum transfer. The experimental PDF 
peak widths as a function of pair distance are extracted 
using the 'real-space' Rietveld program PDFFIT.I13 For 
detailed procedures about modelling PDF spectrum and 
extracting PDF peak widths refer to Ref. 6. 



IV. MEAN-SQUARE RELATIVE 
DISPLACEMENTS IN CRYSTALS 



III. EXPERIMENTS AND ANALYSIS 

The experimental PDFs discussed here were measured 
using pulsed neutrons and synchrotron x-ray radiation. 
The neutron measurements were carried out at the NPD 
diffractometer at the Manual Lujan, jr.. Neutron Scat- 
tering Center (LANSCE) at Los Alamos and the x-ray 
experiments at beam line A2 at CHESS (Cornell). Pow- 
der samples of Ni and a polycrystalline Ce rod were 
loaded into a vanadium can for the neutron measure- 
ments, carried out at room temperature. Powdered GaAs 
was placed between thin foils of kapton tapes for the x- 
ray measurements, measured at 10 K using 60 KeV (A 
= 0.206 A) x-rays. Due to the higher x-ray energy at 
CHESS and relatively low absorption coefficient of GaAs, 
symmetric transmission geometry was used. r- ir- 1 

Both the neutron and x-ray data were correctediJllj for 
experimental effects and normalized to obtain the_t;otal 
scattering function S(Q), using programs PDFgetNllj and 



The PDF peak of simple materials can be welLaxiprox- 
imated by a Gaussian function with a width .E£lEEl The 
mean-square relative displacement of atom pairs, pro- 
jected onto the vector joining the atom pairs, is given 
by 

4 = ([(Ui-Uj).fy]2), (1) 

where Ui, uj are thermal displacements of atoms i and j 
from their average positions .t3E3 The vector fy is a unit 
vector parallel to the vector connecting atoms and 
the angular brackets indicate an ensemble average. This 
equation can be rearranged as 

4 ^ ([ui • fij]2) + ([uj • fij]') -2((ui • fij)(uj • fij)). (2) 

Here the first two terms correspond to mean-square ther- 
mal displacement of atoms i and j. The third term is a 
displacement correlation function, which carries informa- 
tion about the motional correlations. For a monatomic 
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FIG. 2: Theoretical mean-square relative displacement a^j of 
7-Ce as a function of pair distance calculated using Eq.^ and 
the BvK model at 300 K. The inset shows a^j below lij < 
20 A. The solid line corresponds to 2(u^), where (u^) is the 
mean-square thermal displacement of 7-Ce. 
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is expressed in terms of the phonons as 
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-cos(k • ry)], 
(3) 

where Ws(k) is a phonon frequency with wave vector k 
in branch s, n(a;s(k)) is the phonon occupation num- 
ber, ek.s is the polarization vector of the k, s phonon 
mode, N is the number of atoms and M is the mass 
of an atom. As an example, we calculate the afj of Ce 
(Fig. H) using Eq. |[ Ce crystallizes in a simple FCC 
structure (space grouR-J" m3m) at room temperature and 
atmospheric pressurcEil In this calculation, the phonon 
frequency (^^(k)) and polarization vector (ek,s) were ob- 
tained by solving the dynamical matrix using up to the 
8th nearest-neighbor (NN) interatomic force parameters. 
The force. parameters of Ce were determined by Stas- 
sis et a/E3 by fitting the phonon dispepion curves us- 
ing the Born von-Karman (BvK) modelH In Fig. the 
horizontal solid line corresponds to 2(w^), where (i?) is 
the mean-square thermal displacement of Ce. Deviations 
from this line are due to motional correlations in the pair 
motion. The inset shows crf^- below < 20 A, where the 
motional correlation is more apparent. Evidently, motion 
of near-neighbor atoms is highly correlated, and this is 
reflected in narrower PDF peak widths. At larger separa- 
tions {rij > 20 A) the afj values asymptotically approach 
the uncorrelated values because the cosine term in Eq. ^ 
averages to zero. 

As shown in Fig. ^, the motional correlation of atom 
pairs varies significantly as a function of pair distance. 
Therefore it is useful to quantify the degree of correlation 
using a dimensionless correlation parameter which can 
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FIG. 3: Correlation parameter (j> of Ce atoms along vari- 
ous crystallographic directions obtained using the BvK model 
calculations at 300 K. 0=0 corresponds to the uncorrelated 
atomic motion and > indicates that atoms move in phase. 



be defined as follows 
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where erf = ((ui • fy)^). It can be seen from Eq. ^ that 
^ = corresponds to completely uncorrelated motion. 
Positive values of (j) describe a situation where the atoms 
move in phase, and thus the resulting value of cr^ is 
smaller than for the uncorrelated case. Using Eq. ^ the 
correlation parameter (j) can be calculated from the PDF 
peak width measurements as 



(5) 



Figure ^ shows the motional correlations along various 
crystallographic directions in Ce. It is clear that the 
correlation parameter varies significantly with crystallo- 
graphic direction. Along the <110> direction, Ce atomic 
motion shows strong correlation. On the other hand, Ce 
atoms along <100> show almost no motional correlation. 
This behavior comes principally from the characteristic 
elastic anisotropy of cubic crystals.E3 Despite the exten- 
sive orientational averaging of the powder measurement, 
this directional information survives in the data. 

The oscillations in af^ shown in Fig. ^ for Ce are gener- 
ally driven by both the interatomic interactions and the 
crystal structure. This is illustrated in Figs. ^ and ||, 
which show the correlation parameter for a variety of 
FCC and BCC materials. The in Figs. | and |[ are 
generated in the same way as for Ce, that is, using the 
BvK force model, with parameters derived from fits ta 
the phonon dispersion curves found in the literature.c3 
One sees that a common oscillatory behavior in the cor- 
relation parameter is found for all of the FCC crystals 
studied. And similar but distinct oscillatory features are 
observed for the BCC crystals, except Nb. This difference 
in the general behavior of FCC and BCC crystals sug- 
gests that atomic geometry plays a role in the motional 
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FIG. 4: Correlation parameter <j> of FCC crystals Ca, Ni, Ce, 
Al, Au at 300 K, obtained using the BvK model calculations, 
a is the lattice parameter of each crystal. 



CD 

-I—' 

0) 

E 

CD 
1— 

CD 

Q. 

g 

_CD 

<D 
i_ 
L_ 

o 
O 



0.5 
0.4 
0.3 
0.2 
0.1 
0.0 



\/\ 



— I — ' — r- 

Fe 5NN 
Ta 7NN 
Na 5NN 
-NbSNN 




1 2 3 4 5 6 7 8 

r(A)/a 

FIG. 5: Correlation parameter (j> of BCC crystals Fe, Ta, Na, 
Nb at 300 K, obtained using the BvK model calculations, a 
is the lattice parameter of each crystal. 



FIG. 6: Frequency spectrum of a^j for Ce at 300 K obtained 
using the BvK model calculations: (a) INN, (b) 2NN, (c) 
lONN, and (d) uncorrelated far-neighbor atom pair. The area 
underneath the solid line in each figure corresponds to the a^j 
of each pair. 



the broadening of the peak widths. Figure |^ shows the 
frequency spectrum of crf^ in Ce for the INN, 2NN and 
lONN peaks, in addition to that of the uncorrelated far- 
neighbor atom pairs. Referring to Fig. |[ for the INN, 
the low frequency {lojlTi < 1 THz) thermal motion con- 
tributes little to the INN PDF peak width broadening. 
Most of the peak broadening comes from mid-to-high fre- 
quency modes (1 < a;/27r < 3 THz), which contribute 
almost equally to the broadening. This suggests that 
INN pair moves more or less in-phase at low frequen- 
cies, and that the pair motion de-phases somewhat as the 
frequency increases. For the 2NN and higher-neighbors 
peaks, however, where the motion of atom pairs gets 
more and more de-phased as the pair distance increases, 
the medium-range frequency modes predominantly con- 
tribute to peak width broadening. Finally, for the far- 
neighbor atom pairs, where the motion is completely un- 
correlated, the frequency spectrum is spread more evenly 
across all of the frequencies. 



correlations as well as the interatomic interaction. How- 
ever, significant differences in the correlation are evident 
among different elements with the same crystal structure 
as well. For example, in FCC crystals cj) varies from 0.37 
to 0.45 for the INN pair and from 0.05 to 0.2 for the 
2NN pair, fn BCC crystals cj) varies from 0.38 to 0.47 
for the INN pair. Thus, and afj do reflect interatomic 
interactions. 

In addition to a dependence on the atom pair distance 
rij , the a'^j shows an explicit dependence on phonon fre- 
quencies Wfe(s). Therefore, it is instructive to consider 
how phonon modes of different frequencies contribute to 



V. CORRELATED THERMAL MOTION: 
DEBYE APPROXIMATION 

Using the BvK model we have shown that near- 
neighbor atomic motions in crystalline materials are 
strongly correlated. In the BvK model calculation, how- 
ever, the force constants must be known in advance. 
In this section, we simplify the result in Eq. |^ using 
some approximations to describe the effects of the lat- 
tice vibrations on the PDF peak widths-^ithout know- 
ing the forcfi-constants. Following DebyeEZi and Beni and 
Platzmann,E3 we make no distinction between longitudi- 
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nal and transverse phonon modes and take a spherical 
average. Then Eq. || reduces to 

1- 



_ / 2h 



n{uj) 



1 — cos(k • ry) 



(6) 



where (• • ■) is the average over the 3iV modes and N is the 
number of atoms. This result is a general expression for 
all materials rSuid is independent of the number of atoms 
per unit cell.Es Using the Debjte approximation, u) — ck, 
we can write Eq. || as foUowscM 
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sin{u!rij / c) 
Lurij/c 



(7) 

where p{lu) = 'iN{iuj'^ /lod^) is the phonon density of 
states, n(w) is the phonon occupation number, c is 
the sound velocity and uj^ — cki) is the Debye cut-off 
frequency. The Debye wavevector is given by fc^i = 
{GiT^N/Vy/'^ where N/V is the number density of the 
crystal. After integrating over uj, we obtain 
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where $i = /p^°^"^a;(e^ — dx, a; is a dimension- 
less integration variable and Qd (=Swd/A:b) is the De- 
bye temperature. Thi s EesuL t is known as the "correlated 
Debye (CD) model" .E3E3e3 Here, the first term corre- 
sponds to the usual uncorrelated mean-square thermal 
displacements (2(m^)) and the second term is the dis- 
placement correlation function (DCF). The CD model 
was first used to explain XAFS peak widths as a func- 
tion of temperature, and provided reasonable fits to the 
INN and 2NN peak widths. However, this model has 
never been tested beyond the 2NN peak. Here we test 
this simple model against full experimental PDF spectra 
and the BvK model calculations. 

In the CD model the DCF shows explicit dependence 
on the atomic pair distance r^ . The first term of the 
DCF comes from the quantum zero-point motion and the 
second term is temperature dependent. Figure ^ shows 
the rij dependence of the first and second terms of the 
DCF of Ni. The first term of the DCF decreases as l/r?- 
with a cosine modulation. The second term is temper- 
ature dependent and shows a l/r^ dependence. At low 
temperatures (T<^ 9d)-, the temperature dependent DCF 
term is much smaller than the zero-point motion term. 
However, as the temperature increases the second term 
becomes dominant. These results show that the motional 
correlation follows a l/rf^- dependence when T^ Oo and 
a l/vij when T> Od. 

We tested the CD model calculation against the ex- 
perimental and BvK PDF peak widths of Ni at 300 K. 
Figure ^ shows selected experimental PDF peak widths 
and the calculated peak widths as a function of pair dis- 
tance, as well as the corresponding phonon density of 
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FIG. 7: Displacement Correlation Function (DCF) of Ni. 
(a) 1st term of DCF of Eq. |. (b) second term of DCF of 
Eq. I at 300 K, 100 K, and 10 K. Solid line is a fit to a 1/r 
dependence. 



states of Ni. The errors in the experimental PDF peak 
widths were estimated from fitting Gaussian functions 
to the PDF peaks. The error in CD model calculation 
was estimated from the error in the experimental thermal 
displacement. For the BvK model calculation, we used 
up to the 4tlul^N force parameters determined by Bir- 
geneau et a/.. Ell We also compared the 4NN BvK model 
calculations with those of a simple INN BvK model and 
found excellent agreement between them. Our experi- 
mental thermal displacement of Ni {{v?) = 0.00535 A^) 
is about 10% larger than that of the BvK model cal- 
culation. As a result, as shown in Fig. ^(a), the BvK 
model peak widths are shifted downward by roughly 5 % 
overall. The origin of this difference between our thermal 
displacement and that of the BvK calculation is not clear 
but the Debye temperature, determined from our thermal 
displacement using the Debye approximation, ^0=385 K, 
compares well to the specific heat measurement, where 
6»£,=375 KH3 For the CD model calculation the Debye 
wavevector ko = 1.756 was obtained from the atomic 
geometry. The only parameter in the CD model, the De- 
bye temperature, was determined from the experimental 
thermal displacement. 

As shown in Fig. ||(b), the Debye model reasonably 
approximates the 'real' density of states in this simple 
element. However, because the Debye temperature and 
corresponding Debye frequency were obtained from the 
thermal displacement, and not from the sound velocity, 
the low-frequency Debye density of states deviates from 
those of the BvK model by roughly 15%. Nevertheless, 
the CD model calculation of the peak widths shows good 
agreement with the experimental PDF peak widths aij 
except for the overall downward shift. 

We now move to the more complex case of Ce, which 



6 



0.11 
0.10 
0.09 
0.08 
0.07 
0.06 

0.4 
0.3 
0.2 
0.1 
0.0 



(a) 

f 



i t i it 



Experiment 
BvK model 
Debye model 





5 10 

r(A) 


15 


20 


■ (b) - 


BvK model 








Debye model 















4 6 
m/271 (THz) 



, — , 

•< 


22 




0.20 - 












0.18 - 








0.16 - 






& 




l-L, 


0.14 - 


Q 










0.12 - 



(a) 



1.2 

1.0 

0.8 

0.6 

0.4 

0.2 

0.0 
0.0 



f • 



Experiment 
BvK model 
Debye model 



10 15 20 

r(A) 



- (b) - 


— BvK model 












I 



0.5 1.0 1.5 2.0 2.5 3.0 
mlln (THz) 



FIG. 8: (a) Comparison of neutron PDF and BvK model 
peak widths with those of the CD model calculation of Ni at 
300 K. Filled squares: experimental PDF peak widths (only 
a few selected peak widths are shown), Filled circles: BvK 
model calculations, Dotted line: CD model calculation using 
Debye temperature ^d=385 K, Debye wavevector fc_D = 1.756 
A~^. (b) Solid line: Ni phonon density of states calculated 
using BvK model. Dotted line: Debye density of states with 
the same area as the BvK calculation. Debye cut-off frequency 
ujD/2-n-=8 THz. 



FIG. 9: (a) Comparison of neutron PDF and BvK model 
peak widths with those of the CD model calculation of Ce at 
300 K. Filled squares: experimental PDF peak widths (only 
a few selected peak widths are shown). Filled circles: BvK 
model calculation. Dotted line: CD model calculation with 
Debye temperature 0d=117 K, Debye wavevector fc_D = 1.1986 
A~^. (b) Solid line: Ce phonon density of state calculated 
using BvK model. Dotted line: corresponding Debye density 
of states, using Debye cut-off frequency ujd /'2'it— 2.44: THz. 



needs long-range forces to explain anomalies in the dis- 
persion curves. Up to the 8th NN interactions are re- 
quixed to reasonably fit the phonon dispersion curves of 
CeE3 Fi gure y shows selected experimental PDF peak 
widths and calculated peak widths as a function of pair 
distance, as well as the phonon density of states of Ce. In 
this case, the BvK calculation of the peak widths shows a 
good agreement with the experimental PDF peak widths. 
We also found that a simple INN BvK model calculation 
shows very good overall agreement with that of 8NN BvK 
model, except for the 2NN PDF peak width which devi- 
ates by ~ 3 %. For the CD model calculation, a Debye 
wavevector ko— 1-1986 and a Debye temperature 

6d = 117 K (at 300 K) were obtained from the atomic 
geometry and thermal displacement of Ce ((w^)=0.0231 

) , respectively. Even for the more complex system like 
Ce, the CD model calculation of the PDF peak widths 
reproduces the overall -dependence rather well, except 
for a few detailed features. 

We now consider the T-dependence of aj, for specific 
atom pairs of Ce. These are shown in Fig. |o| for the first 
few NN pairs, as well as the uncorrelated atom pair. All 
of these curves exhibit characteristic Debye-like behav- 
ior, i.e., linearity in T at high temperatures, but curving 
over to a common zero-point value at T = 0. In general, 
the CD model calculation of the temperature dependence 
of the erf - shows good agreement with that of the BvK 



model calculation, except for the 2NN pair (Fig. |l^). Re- 
ferring to Fig. ^ we see that BvK calculations of the afj 
of the INN and 3NN lie very close to the Debye predic- 
tion, but the 2NN is displaced significantly upwards. The 
deviations between the BvK and Debye models increase 
as the temperature increases. Apparently, if the af^ lies 
on the Debye prediction at one temperature, the Debye 
model will also predict its temperature dependence cor- 
rectly. Conversely, the temperature dependence of afj 
will be underestimated or overestimated depending on 
whether it is displaced above or below the Debye curve 
at the lowest temperatures, respectively. This is at least 
true for Ce. 

As a final example, we compare the CD model calcula- 
tion of PDF peak widths with those of GaAs determined 
experimentally. This is distinct from the above exam- 
ples due to the presence both of more than one atomic 
species and directional covalent bonding. We also com- 
pare the CD model calculation with a lattice-dynamics 
calculation using the Kirkwood potential, which has been 
shown to be a good b asis .for describing semiconductor 
compounds and alloys .liaE^ Figure |ri](a) shows the atom 
pair dependence of the PDF peak widths of GaAs at 
10 K. In the Kirkwood model, the potential parameters 
are obtainedJajtitting the nearest and far neighbor PDF 
peak widths .EjS In Fig. |ll|(a), the Kirkwood model cal- 
culation using bond stretching (a=96 N/m) and bend- 
ing (/3=10 N/m) force constants shows quite good agree- 
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FIG. 10; Mean-square relative displacement a^j of Ce as 
a function of temperature. Upper panel: nearest neighbor 
(NN), second NN (2NN). Lower panel: third NN (3NN), fifth 
NN (5NN) and DW. DW represents uncorrelated far-neighbor 
pair. Symbols are the Born von-Karman (BvK) model calcu- 
lations and lines are the corresponding CD model calculations. 
In BvK calculation, the Debye temperature, 6'£) = 117 K is de- 
termined at 300 K. 



ment with the experimental PDF peak widths. In the 
CD model, the average mass of Ga and As is used. The 
Debye wavevector kn = 1.382 and Debye temper- 
atm'e 6*15=250 K were obtained from the atomic geom- 
etry and by fitting the far-neighbor PDF peak widths. 
In the CD model calculation, the near-neighbor peaks 
(r < 5 A) are ~ 10% broader than those of the ex- 
perimental peaks. Referring to Fig. |l^(b), the Debye 
model does a poor job of describing the GaAs phonon 
density of states; for example, the high frequency optic 
modes, 6 < a;/27r <8 THz, are totally missed. Instead, 
the Debye model over-estimates phonon modes between 
3.5 THz < w/27r < lojj. This poor description of the 
phonon density of states, as well as the use of the average 
mass of Ga and As for M in Eq. |^, leads to deviations in 
near-neighbor peak widths from those of the experimen- 
tal peaks and the BvK model calculations. Therefore, 
the deviations in CD model calculations basically reflect 
the limitation of the CD model in describing motional 
correlations in a system with optic modes. Nevertheless, 
the CD model, with a single parameter 9d determined 
from the thermal displacement, serves as a good first or- 
der approximation to the PDF peak widths, even in more 
complex systems like GaAs. 



VI. DISCUSSION 



The mean-square relative displacements afj and the 
corresponding correlation parameter shown in Figs. 0, 0, 
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FIG. 11: (a) PDF peak widths of GaAs as a function of 
atom separation at 10 K. Filled squares: Experimental x-ray 
PDF peak widths, Filled circles: lattice dynamics calculation 
using the Kirkwood model (a=96 N/m, /3=10 N/m), Dot- 
ted line: CD model calculation using the Debye temperature 
00=250 K, Debye wavevector fcD = 1.382 A"^ (b) Symbols: 
GaAs phonon density of states calculated-iusing local-density 
approximation density functional theoryu. Solid line: Debye 
density of states with the same area as the LDA calculation. 
Debye cut-off frequency lo-aj^-n — 5.22 THz. 



^ and H, H present two interesting pieces of information 
about the atomic motions in a crystalline material. First 
of all, they show that nearest-neighbor atomic motion is 
significantly correlated. Second, the details of the mo- 
tional correlations as a function of pair distance display 
structures which deviate from the predictions of the sim- 
ple CD model. Here we can raise some interesting ques- 
tions. How is this structure in the motional correlation 
of atom pairs related to the underlying interatomic po- 
tentials? Can one extract the potential parameters using 
an inverse approach to model the PDF peak widths with 
the potential parameters as input? 

Reichardt and PintschoviuJl argued that the calcu- 
lated PDF peak widths as a function of pair distance 
are rather insensitive to the details of the lattice dynam- 
ics models used to calculate afj. They found that PDFs 
calculated using either very simple or complex models 
didn't show significant differeii|Ces. A similar conclusion 
has been reached by Graf et al.,i3 in contradiction to pre- 
vious claims by Dimitrov et al.n Indeed, the magnitude 
of errors implicit in the measurement and data analy- 
sis appear to be comparable to the effects that must be 
measured to obtain quantitatively accurate potential in- 
formation using this approach.EI The conclusions of Re- 
ichardt and Pintschovius and Graf et al. and Thorpe et 
al. are largely borne out by the present work; e.g., the 
grossly oversimplified CD model, which neglects elastic 
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anisotropy and parameterizes the dynamics with a single 
number ^d, is rather successful at explaining the smooth 
Tij-dependence of the PDF peak widths. 

Thus, when the BvK force parameters are not avail- 
able, we have shown that the correlated Debye (CD) 
model is a reasonable approximation to describe both the 
smooth r.y -dependence and the temperature dependence 
of afj in simple elements. Considering the poor corre- 
spondence between the Debye phonon density of states 
and the BvK density of states, the reasonable agreement 
between the BvK model calculations of tr? and that of 
the CD model is rather surprising. This confirms that 
the PDF peak width is rather insensitive to the details 
of the phonon density of states and the phonon disper- 
sion curves, as suggested by Reichardt and Pintschovius 
and by Graf et ai. Any information about the inter- 
atomic forces in the PDF peak widths is contained in the 
small deviations of the cr?- from those of the CD model 
calculations. Therefore, extracting interatomic potential 
information from the PDF peak widths is unlikely. How- 
ever, these deviations could possibly yield some average 
phonon information. For example, recent calculations by 
Graf et alx^ show that one can obtain phonon moments 
within a few percent accuracy for most FCC and BCC 
crystals using the nearest-neighbor force parameters ex- 
tracted from a theoretical BvK PDF spectrum. This re- 
sult indicates that the PDF spectrum contains some av- 
erage phonon information, although it doesn't provide 
detailed phonon dispersion information. The average 
phonon information, such as phonon moments from the 
PDF peak widths, will be a useful complement to optical 
and acoustic techniques that yield zone-center informa- 
tion in situations where single crystal measurements are 
not possible. This complementarity also extends to the 
extraction of Debye- Waller factors from powder diffrac- 
tion measurements. 

Finally, a comparison of the CD model calculations of 
the PDF peak widths in GaAs with those of experimental 
PDF and Kirkwood model calculations shows additional 
limitations of the CD model. In the CD model calcula- 
tion, the near-neighbor PDF peaks below r < 5 A are 
about 5-10% broader than those of experimental PDF 
peaks. This is due to the poor description of GaAs 
phonon density of states by the Debye model. Since 
the sine term in Eq. |^ over- and under-weighs certain 
phonon modes depending on their frequencies, the re- 
distribution of GaAs phonon density of states in a realis- 
tic model causes deviations in near-neighbor peak widths 
from those of the CD model. One way to improve the 
model calculation in materials which have optic modes 
might be a hybrid model that combines the correlated 
Debye and Einstein models. Such a—hybrid model has 
worked quite well in the case of AglEj 
VII. SUMMARY 

In this paper the mean-square relative displacements 
{erf A of atomic pair motion in crystals have been stud- 



ied as a function of pair distance and temperature using 
the atomic pair distribution function (PDF). The experi- 
mental PDF peak width and the BvK model calculations 
of afj as a function of pair distance show that the near- 
neighbor atomic motions are strongly correlated. The 
extent of these correlations depends both on the inter- 
atomic interactions and crystal structure. Thus, a proper 
accounting of the lattice vibrational effects on the PDF 
peak widths is important in order to better understand 
the effects of static and dynamic disorder on the PDF 
peak widths in disordered systems. Details of the PDF 
peak widths vs. seen in the BvK calculations are 
well reproduced in the measured data indicating the ac- 
curacy of the measurements. Most of these details orig- 
inate from the elastic anisotropy of the crystal which is 
especially apparent in FCC crystals. We showed that the 
CD model reproduces the average features of the lattice 
vibrational effects on the PDF peak width with just one 
parameter, which is determined from the measured ther- 
mal displacement (m^). Therefore, this simple model can 
be used as an important adjunct when using PDF to ex- 
tract static and dynamic disorder disorder information 
from materials with local lattice distortion. In addition, 
the T-dependence of the CD model largely agrees with 
the BvK model calculations. Good agreement between 
CD model and experimental PDF peak widths indicates 
that the PDF peak widths are rather insensitive to the 
details of phonon density of states and the phonon dis- 
persion curves. Any information about the interatomic 
forces in the PDF peak widths is contained in the small 
deviations (< 5%) of the cTy from those of the CD model 
calculation. This makes the extraction of interatomic po- 
tential information from PDF peak widths alone unlikely. 
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